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ABSTRACT 


The propagation of simple acoustic transients in an 
isovelocity water layer was investigated for the purpose of 
studying a relationship between the propagation of transients 
in ducts and in layers. The Laplace transform method was 
used to obtain approximate solutions for the acoustic 
pressure, particle velocity, and particle displacement result- 
ing from a step-function input in velocity and the particle 
velocity resulting from a gated sine-wave input in velocity. 
Computer programs were written to evaluate and graph 
the resulting waveforms. The waveforms resulting from 
step~function and gated sine-wave inputs in velocity were 
observed and compared with the predicted waveforms. With 
the use of Mylar transducers, good correlation was obtained 
in the rege of validity of the theoretical solutions. It 
was found that previously published solutions for the 
propagation of transients in ducts became the leading terms 


of the associated solutions for the layer. 
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liye INGER ODE cCTiIen 

While most investigations of the propagation of energy 
in waveguides have been done for electromagnetic cases, 
some have considered the propagation of pulses of acoustical 
energy in waveguides! 122:3,4) , The dispersion of an ultra- 
sonic pulse propagated through a waveguide of circular cross- 
section filled with water has been used for the approximate 
determination of the frequency spectrum of the original 


(1). 


pulse It has been show. that the dispersion of a 
rectangular pulse along an acoustic waveguide of rectangular 
cross=Section can be used to calibrate the waveguide so it 
can be used to predict the effects of similar dispersion on 
the transmission of pulses of other shapes that are not 


(2). 


amenable to actual computation Pulse=compression 
studies in an acoustic waveguide show that a frequency 
modulated input pulse can be compressed into an output pulse 
of shorter duration and of higher peak-amplitude!3). in 
addition to the waveguide experiments, investigations of 
transients in layered media have been made! 226) . Wie. 
peruse” investigated the propagation of a pressure 
transient in a surface layer of liquid which lay over a deep 
bottom layer of liquid. Many of these studies have used 
approximation techniques such as the method of stationary 
phase. A summary of these techniques useful for calculating 


the resultant propagated waveform has been presented by 


iS, 


J. R. Wait!7). The original work of Rayleigh on scatter- 
ing and radiation of harmonic disturbances by small 
cylinders was extended to transient disturbances by J. W. 
Miles using the Laplace transform'8), R. B. Barakat 
used the Laplace transform to investigate the radiation of 
acoustic pulses of various shapes from an infinitely long 
circular cylinder‘ ?? - 

Recently, exact solutions for the propagation of simple 
transients in waveguides (step-function and gated sine-wave 


(4,10,114). Some work has also 


inputs) have been published 
been done which suggests a correlation between these exact 
solutions in ducts and solutions for the propagation of 
acoustic transients in isovelocity laos 

The purpose of this research is to investigate the 
region of validity of the suggested correlation and to see 
what, if any, restrictions must be placed on the resulting 
solutions describing the transient response of the layer. To 
facilitate this comparison, those results of the propagation 


(4,12) 


of transients in ducts which will be of value are 


summarized in Appendix A. 
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ae AR wa Gise 4 
Development of Solutions 
The response of an isovelocity layer to simple transient 
excitations will be developed using the following assumptions: 
(a) Cylindrical velocity=source of radius a. 
(b) Constant layer depth, d. 


(c) Perfectly-Reflecting upper and lower boundaries. 


2 





An Isovelocity Layer 


Figure 2.1 
For radially-symmetric propagation in the layer indicated 
in Figure 2.1, the wave equation for the acoustic velocity 


potential b can be written in cylindrical coordinates as 
(= + to + oO _ ! a* O( + O (2.1) 
Bre? or Or >> ce pur Se ) = 


where r is the radial direction and z is the transverse 


direction. 


LS 


The separation of variables 
b airit)= dy, pst) = 2, 2,3) Bont) (2.2) 


leads to two partial differencial equations: one for a 
“transverse” function Ee which depends only on the height 


z of the layer 
i 
Od, | 
[a ae = eae 
[33> (FS) 122g) = O , ins 
and one for a "radial" portion 9, of the acoustic velocity 
potential, 


z n\~ a 
Lda td -(S)- a ge[ainn-O | 


a 


which depends on both time, t, and the distance from the 
origin, r. ‘The value of the separation constant @, depends 
on the nature of the boundaries at z = 0 and d. For 
example, if we have a presSsure-release top and bottom, 
then p --eof =O at z = 0 and d and Eq. @&.3 has the 
Ot 

solution 

Z,,(3) = Dsip On ( 2 

c 

where Gy, = nTrc/ 4 . Furthermore, since the source is 
assumed symmetric with respect to the plane z = d/2, just 


those modes having pressure antinodes at da/2 will be excited. 
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Thus only odd values of n are allowed. ‘Therefore Eq. 


2.5 can be rewritten as 
Z,,(3) = Den (27 2) ae 1,3, 5, ete, (2.6) 


Since this transverse function affects the form of Have. 
only through the allowed values of O) it is convenient to 

establish a set of reduced acoustic variables which suppresses 
the function z,,(3) 


The well-known acoustic relations 


7, > = aon 5 EL) 


pe 
Lu, = VS (2.7b) 


and 


= 2. o. Cea) 


where n is the mode number, Py the acoustic pressure for 
the nth mode, 3k, the particle velocity for the nth mode, 
and Bn the radial particle velocity for the nth mode can be 


msed Withe@q. -£.k tousmow, that 


Pr = -¢ 2ng) 2 (nt), om.) 


dle, 


The "fragment !"' ne of the acoustic pressure for the nth 


mode is defined as 


— tr - 
ie — M.,, a Tat _ (2.9) 


where Z. = 2. (r,t) is the radial portion of the velocity 
potential for the nth mode which was previously defined by 
Fiq. 2-2. Combination of Eas.2.2 and 2.7c allows the 
definition of the "fragment" of the radial particle-velocity 


for the nth mode, 


— LL my, — } - 
aie — Janie ‘Sa ©, (2. i0% 


Since the transverse function a (z) is independent of 
time, the radial particle-displacement for the nth mode may 


be obtained by integration of u,, with respect to time. This 


n 
enables the "fragment" of the radial particle-displacement for 


the nth mode to be written as 
| € re 
§ = n Sey ip, dé (2.11) 


These reduced acoustic variables or "fragments" $,, Ur py: 
Pa» and ae will be used throughout the remainder of this 


development. 


iS 


Since we are considering the propagation of transients, 
Eiq. 2h, a differential equation in two unknowns, cannot be 
solved by separation of variables. The Laplace transform 
method will be used to obtain the desired solutions. 

The Laplace transform of Eq. 2.4% with quiescent 
initial conditions results in a differential equation in one 


unknown which may be written as 


iS 2. 2 — 


where the bar over the function denotes the Laplace trans-— 


form of itmhnageeimction : 


B (r,s) = L fg, cna} : f ot D (r,t)dt, 


The change of variable 


Gee. Vee? (2a 
C 


gases Ea. c.le into the form 


a Z. | > g — 
<= ‘Gage =O a 
Oe Je Of. 

The solutions to Ea. 2.8 are I, (pe ) and K VA) ila 


modified Bessel functions of the first and third repels 2) 


Conservation of energy requires D,. to go to zero as Sx~@, 


Ee 


so that the only acceptable solution of this differential 
equation is K, (fr). The Laplace transform of @, can 


now be expressed as 


Z = C Kip) (2.15) 


where C! is an amplitude constant to be determined from 
the boundary conditions. 
With the use of the relationship expressed by Eq. 2.7c, 


the Laplace transform of eo becomes 


vn = =C!2Ko(p) 30 i. mr (2.16) 
ae OA. sce =—C Ist = 





The known boundary condition at the source, located at 


r =a, for propagation in the nth normal mode, 


Ven (o,t) 7 oF, (2,2) = V(t) (ee 
or ) 


has the transform 


ny. (at) = Va. (2.18) 
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Equations 2.16 and 2.18 allow evaluation of C!, 


/ — Voc 


1 rr K (S\seF") ) (2.19) 





so that 74 rn Can be expressed as 


Vr = Vo. K, (E\s*+on” ) ° 
K, je [s*+o.~ ) 


(2.20) 





The substitution for Clin Eq. 2.15 allows J, to be written 
as 


= ai. K.(ZVs%on) (22m) 


Po® ewer \ (Von) 


The expression for the Laplace transform of the fragment 


of acoustic pressure, p Can benoUtaimecd trom iq... cat 


rn’ 


if we recall that Aer = -¢ 286 corresponds to Pn =-e5, 


in the transform domain. Thus, 


HY, = ecsV. k (EVs) (mee) 
rn Ree a TF 
(sor K(e¥s¥o7) 


ra A 


The Laplace transform of the particle-displacement 
mira cm ent, 8 rn » can be obtained from VU" nn by either 


integration with respect to t in the time domain or multi- 


plication by 1/s in the transform domain. The result is: 


Se. ie 
es = VS 4? 
ee = Vo KEYS or) e (2.238 


Inverse transforms of Eq. 2.20 — 2.23 have not been 
tabulated, so approximation methods must be used. 


dao KA. (Eq. 2.13) is greater than unity the asymptotic 





form 
a dy? = 
Ky (4) = =a eo (| + 5 +--+) (2.2h) 


ele, 13))m 


may be used The requirement AA > 1 can be 


satisfied by restricting c/a), to values less than unity. 


-UbStitution of Esq. 2.2h intomiiGd.s2.20 yields 


ai \ wT > 
Ue, = Vv. 2/1 etea,h exe “Esme |(I+gpuar) (2.25) 
3 
exe 2/55] ((+ Fa) 








ee 


which simplifies to 


Ven ~ Ve exe|- TY stun )[ +30 rain 

















MB eley) 
[\* caper] 
Qa ; 
Cie $43; 
where T = r-a is the time-of-flight of the signal from 
e 
source to receiver. Since c/aw, is restricted to values 
less than unity, the denominator can be expressed as 
| Oy 5 
— =a a 
| ; 2 | BEVsuw2 ) O25 


so that we have 


Ven Teen VAS exe[-TW2r] [ [i+ § sty | |1- ge ecm | (2.28) 


which simplifies to 


(eo .=\ exe|-T\ S40, | -\, 3c ee (eS exo- ( TVs? | ( 222g) 


8 BoNsau> Wrr 


Vd 
where terms of order (c/a w) have been dropped. 
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With similar manipulations the other desired acoustic 


variables can be obtained from Eqs. 2.21 = 2.23: 


(EG,- he a ee Be Sch (in 8) xr Er Rear] (2.30) 


Venue W,> 


(E %= es Vy exelWsror] — sects Ve oV. ((+ #) exe (-TYs Oe | (2.31) 
(seo> Ba(s¥w?) 





(2.32) 
Bo\istw,> 


fe € = Ve. exe(-Tsor | _ Ye (-) exef-T VStar | 


Equations 2.29 = 2.32 approximate the Laplace transforms 
of the desired acoustic variables for the boundary condition 
Va(t) at r = a with the restriction c/aw - Za, 

Comparison of Eq. 2.29 with Eqs. A.7/ and A.9 shows 
that the terms of Eq. 2.29 have known inversions. F'or a 
step input, for which V, = “7 the inverse transform for 
the fragment of the radial particle-velocity can be written 
down directly from the known results of the duct (see 


Appendix A 


; 2) a 
(E Va= J cone) + 2B AE) J toe) — 2 (1-2) 22 ($5 * J nes) (2.33) 


au 


where ARG =e \ ee and it is always understood that 


all Se) icLeme Vanish fon eeu 

The inverse transforms of the other acoustic quantities 
can be found similarly by using the known transform pairs 
from the duct case and the additional transform pairs 


developed in Appendix A, ‘The results are 


mi 2 nenee 
le 6° 232)" J (reo — 2E( raya (SB) J | (AR9 (2.34) 


A(2m-) 


& I 
[op = ¢¢{J,faae) ~ Se -( ie DO yer) * J (nee) (25355) 


£6 2] omt) Toe) 38 (yeah Jer | 0 


The quantity t' = t - T appearing in Eqs. 2.33 - 2.36 and 
in ARG as (t = T)(t + 7) = fe _ ~* represents the time 
which has elapsed since the signal traveling with velocity c 
first arrived at the position r. It will be referred to as 


the delayed time. 


Equations 2.33 -— 2.36 are the solutions for the transient 


a5 


response of an isovelocity layer to a step input in velocity 
atr=aandt = 0. 

Similar procedures may be used to obtain the fragment 
of the radial particle-velocity for a gated sine=-wave input 


See 


in velocity for which = ezpust . 


lc Ven U bowe,nee) +U (ws. ARG) + 24( +8) Updo) (253 


where Uru, nee) = Lommel function of two areuments\—/ eee 


A Virm 
Uy Grs,ned) = Zoe (= an .) 


and mrss W(t-TY/ 4S), 


Other acoustic quantities could not be completely solved in 
this case. Details are shown in Section A=-3 of Appendix A. 

Comparison of the above solutions with the solutions 
for the duct listed in Appendix A shows that in all cases 
the solutions previously obtained for the duct become the 
leading terms in the solutions developed for the layer. 

A physical understanding of this can be gained by 
comparing propagation processes in a duct with those in a 


layer. While a duct presents a constant cross-section to 
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the propagating wave form for all time, a transient 
propagating outward in a layer is presented with an 
effective cross-sectional area which increases as re How- 
ever, for sufficiently large r the cross-sectional area will 


not change appreciably as the waveform propagates through 


unit distance. In the limit, therefore, we expect the 
radial divergence contained by the term = x in Eq) '2e5 


Be. become insignificant. In this limit, propagation processes 
can be described by an equation virtually the same as Eq. A3 
which describes propagation in a duct. 

For this reason it is not surprising that the solutions 
for the duct become the leading terms of the layer solutions. 
The remaining terms in the layer ecole are correction 
terms necessitated by the fact that the layer does not 
actually present a constant axonensecaen to the propagating 
wave form. 

The solutions for the transient response of an 
isovelocity layer have been developed with the restriction 
c/aw n< i: Since the cutoff frequency is given by 
WwW, =n Cf » this restriction places a limit on the 
relation between the depth of the layer and the size of the 
source, d/a nit . If it is desired to remove this restric-— 


tion a different expansion of the K., (n) 's involved must be 


used. 


a 


af nvestigation of Solutions 


The solutions developed for the layer have been obtained 
for a single mode. In determining the contributions to the 
solutions from other modes, attention must be given to 
the transverse behavior of the transmitted and received 
signals. ‘The boundary condition on velocity at the source 
at r =a says that V,(t) vanishes for z €0 and z >d and 
is constant for 0<€ z<¢d. This condition can be met by 
representing Vif t) by a spatial Fourier-superposition of 
| u(z) = Bil z—d) | in terms of the normal modes across the 


face of the source: 


| 


Notice that the higher modes contribute as oa for this 
input. ‘These higher modes have positive and negative 
portions in their waveforms over the depth of the layer. 
These positive and negative portions will partially cancel 
over the face of a receiver which extends over the layer 
depth. In effect the receiver will average the received 


(1) 


Signal over the area of its face 


d Pa os 
a { ( ) (2n-1)i73 (2.39) 
i TT i} 2 


n= 
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Evaluation of the above integral gives 


2n-\ 


n= bd F 


rw a 
mi 


[sh (280) 


Therefore each mode contributes as (xt) or 1, 14eF 
1/25, 1/49, for the first four excited modes. 

A computer program was written to evaluate the 
effects of these higher modes on the solutions developed. 
Computer predictions considering (a) the fundamental mode 
plone and (b) the fundamental mode plus the next four non= 
zero higher modes, were plotted versus delayed—-time. 
Comparison of these two predictions as presented graphically 
by the computer showed changes in the amplitude of less 
than 1 per cent and shifts in axis crossings of less than 
0.01 local wavelengths. Since these effects are negligible 
compared to the experimental errors encountered in this 
research, it was sufficient to retain only the fundamental 
mode (n = 1). 

Examination of Eq. 2.36 shows that it consists of 
momeerms. Lhe first (main) termi1s from the duce 
solution and the second (correction) term is from the layer 
effect. An estimate of the effect of the correction term 


can be made by combining the leading terms of the indicated 


Bo 


sums. Equation 2.36 can be rewritten as 


4 


e £-T\/2 
ai aa J (are) = (-2( 5) J, (are) (Zale) 


Under the assumption t' = t-T<€ T, which limits considera- 
tion to the leading portion of the transient waveform, we 


have t&T, and Eq. 2.41 becomes 


(= YY ay ne : 
a 2 = = Je) Saisie) aR (2.02) 
where 9 = ARG = wo, \lt'2T 


VY 
and Bt) = Fos,(2)G5) * 


Substitution of trigonometric approximations of the Bessel 


(13) 


functions in Eq. 2.he gives 


Teer WTR cos (©- 37) - BE) acs(6-S77) (2.13) 


By combining the cosine terms of Eq. 2.43, we find the 
expression for the fragment of radial particle displacement to 


be 


Bear wT = V+ acer cos (6-31 +) (2.40) 


30 


I(t). 


where “ = tan 

Equation 2.44 shows that the correction term in the 
layer solution introduces a phase angle x that was not 
present in the duct solution. ‘The zeroes of the right hand 
side of Eq. 2.4) give the times for which the particle 
displacement vanishes. It has been shown that the duct 
solution becomes the leading term of the layer solution; 
therefore, the zeroes of the first term of Eq. 2.43 give the 
nulls, or axis crossings, of the particle displacement in a 
duct. If these times are evaluated, the shift in time of 
associated ant crossings between the layer and the duct 
solutions can be predicted. 


The zeroes of the layer solution, designated by (a 5 


are obtained when 


cos(6-3%+¥) =O, (2.45) 


or 


S- a +X — (222!) Ty mnahasete, (2.46) 


Sis 


The zeroes of the duct solution will be designated by ee 


They are obtained when 


cos (©- 377) =, (2.47) 


or when 


(o-2) =. (Heat ir m =42,3, ek. (2.48) 


Since we have t'<¢ T, c/aw ,< 1, and r= a, so that 


Be ie 
ton ¥ = BE) = 25, (1-2)\ A CC 


a small-angle approximation can be apphed with the result 


{ 
ton Y= ¥=eTT 2a 


where 


oye 


Substituation of this expression for x ino. - qumecomand 


evaluation of © in Eqs. 2.46 and 2.48 gives 
(9, Var } \t, 2 = (= are (2.51) 


zisalte! 


/ — 
(wT) ee = (2+)? (2.52) 


Manipulation of Eqs. 2.51 and 2.52 reveals that the ratio 


t 


“1/1. can be expressed as 


* 1 
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predicting the shift in axis crossings. 

This result is approximate, since only the first terms 
of the sums in Eq. 2.36 were used to develop it. To 
evaluate the accuracy of this approximate result, computer 
programs were designed to evaluate the predicted waveforms 
of Eqs. 2.33 and 2.35 = 2.37. They were written so that 
the main terms and aS correction terms of these equations 


could be evaluated separately. ‘This permitted comparison of 


2), 


the layer solutions and the duct solutions and showed the 
behavior of the correction term as a function of c/aW 7 
These programs, one for the step input (Appendix A) and 
one for the gated sine-wave input (Appendix D) consider 
up to 11 terms of the sums involved. ‘The programs give 
printed and graphical output of the main terms, the 
correction terms, and the complete expressions of the 
equations mentioned above. Figure 2.2 shows the velocity 
Pet weoMtPiodlicted by Che cOmpULCm umm E 1oMreo ue .] morn 
show the complete expression, the main term, and the 
correction term of the displacement fragment for several 
values of c/aw a 

Information about the envelope of the displacement 
fragment waveform can also be obtained from Eq. 2.44 


\ a 
which gives the amplitude as \+ 8 (t) . Under the 


assumption t! <€¢ T, B(t) becomes 


4 ¢ ae (2.5) 
aey= 25, (SVS 


and 


/ 


R(t)” = = (a, (Sy To (215e 
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Therefore, the amplitude of the initial portion of the 


waveform under consideration is 





The approximate solutions for the propagation of simple 
transients in an isovelocity layer have been obtained, and 
computer programs to evaluate these solutions and predict 
the axis—crossings and amplitude=modulation in the layer 


meave been written. 
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(b) Layer correction term 


Fracment or the radiel particle displacement, Sons ves delayed time, 
t', as predicted by the computer program. c/ao,= 1.00, r= 61.5 cn, 
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Figure 2.4 
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3. EXPERIMENTAL APPARATUS 

‘The experimental apparatus used in this research 
consisted of a tank, several transmitters and receivers, 
and the associated electronics equipment. ‘The laboratory 
Setup |iS@shewm in Fic wean, 
Tank 

‘The water layer was contained in a tank constructed 
of 3/4 inch marine plywood which was 30 cm high, 1 m wide, 
and 1.5 m long. ‘The plywood was coated with Varathene 
to provide waterproofing, and lined with 1/) inch 
polyethylene foam. ‘This prowided the desired pressure 
release bottom boundary for the layer. The water-air 
interface provided the top pressure release surface. With 
the sides and ends of the tank also lined with polyethylene, 
the tank could be excited as a pressure-release cavity which 
served to accurately determine the depth of the layer. 
This technique is described in more detail in section lh. 

Cylindrical sources could not be used because the tank 
was too small. ‘Transmitters were therefore built in the 
shape of semi-cylinders and mounted on a rigid plate. A 
neoprene-coated steel plate, 3/8 inch thick, was made to 
fit in one end of the tank to act as the rigid boundary. 
This rigid surface served as an Bee a) Mirror.) SOmien aie 
the reflected pressure appeared to be emitted from an 
image source (an identical semi-cylinder back-to-back with 
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the real source). Under this condition cylindrical symmetry 
was approximated. Figure 3.@é illustrates the geometry 

of the layer. 

‘Transducers 

Mylar dielectric transducers were selected to be used 
in the layer for the following reasons: 

(a) Their broad-band transmitting and receiving 
characteristics allowed the transmission aie reception of 
the transient excitations without introducing significent 
frequency=—sensitive modulation of the signals. 

(b) They could be conveniently constructed to the 
required semi-cylindrical shapes and nee. 

(c) Previous Wonk eae has indicated that, when 
driven with a low-impedance voltage source, they generate 
an acoustic signal for which the normal particle velocity at 
the source is proportional to the input voltage. 

Two transmitters were constructed to the required 
semi-cylindrical shape, using the techniques described go 
one with radius 1.0 cm and one with radius 2.5 cm. Both 
transmitters were 12 cm high which corresponded to the 
maximum layer depth to be encountered. These dimensions 
were chosen in order to cover the desired range of values 
for c/a ,. 
A Mylar dielectric receiver, identical in nature to the 


transmitters, but with a plane, rectangular face measuring 
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3 cm x 14 cm, was used to measure the averaged (over 


depth) radical particle velocity and displacement in the 
meee By proper adjustment of its electrical load, the 
receiver could be made sensitive to saver particle velocity 
or displacement. When terminated by the 1 megohm 
impedance of the voltage amplifier, it was displacement 
Sensitinee when shunted with an impedance of 100 ohms, it 
became sensitive to velocity. ‘Transmitters and receivers 
are shown in FPiigure 3.3. 

Two pressure-Sensitive receivers were constructed, 
but neither proved to be satisfactory. A probe using a 
1/8 inch barium titanate cylinder as the aenziac element 
was built, but its low sensitivity resulted in an unacceptable 
signal to noise ratio. A 2.5 cm x 2.5 cm sensing element 
cut from a 0.5 cm Glennite disc was also constructed. 
However, the resonant frequency of the element was 
SUP Picieaelyy close to the frequencies at the front of the 
received waveform so that ringing of the element caused 
severe distortion of the received signal. For this reason, 
and because it was felt that the velocity and displacement 
measurements would adequately describe the waveforms, 
pressure measurements were not taken. 
Bilectronies Equipment 

Figure 3.4 shows the block diagram of the equipment; 


Table 3.1 identifies the items in the figure. ‘The oscillator 
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Table 3.1 


Electronics Equipment 


OSC Audio oscillator, TS-382 D/U 

TBG Tone burst generator, General Radio type 1396-A 
Oe Unit pulse generator, General Radio type 1217-A 
ae =o Filectronic counter, Hewlett Packard 5233-L 


PWR AMPPower amplifier, Hewlett Packard )67-A 
AMP i Voltage amplifier, Teletronix type ileal 
AMP 2 Wide band chain amplifier, SKDL model 202-D 


RC FIL Homemade variable low-pass RC filter with 
400 kH, roll-over. 


SKL FIL Electronic band pass filter SKL model 302 
SCOPE Dual beam oscilloscope, Fairchild model 766-H 


LBN Load and bias network (made for earlier work) 
Shown in Pig. 3.5 


7 


and tone burst generator were used as the initial signal 
source for the gated sine-wave input. For the step input, 
the unit pulse generator was used. 

The received signal was sent from the receiver 
through the load and bias network, amplifiers, and filters 
to trace "B" of the dual beam oscilloscope. The input 
signal taken from the output of the power amplifier was 
displayed on trace "At!, Measurements of the received 
waveform axis crossings and amplitudes were made directly 
from the oscilloscope display, or photographed for later 
analysis. Triggering of the scope was provided by gating 
signals from either the tone burst generator or the unit 
pulse generator. ‘The load and bias network shown in 
Figure 3.5 served to supply the. 300 volt bias needed by 
the Mylar transducers without allowing this voltage to 
appear elSewhere in the circuit, and to provide the proper 
electrical load for the receiver so that either particle 
displacement or velocity could be observed. The 22 ohm 
resistor in the transmitter circuit was needed to suppress 
parasitic oscillations caused by the transient electrical 
excitation of the capacitance of the transmitter and the 
residual inductance of the power amplifier. 

The rise time of the system was experimentally 


determined to be about 0.5 Ad Sec. 
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lee Bete Pe LMeN TAIL PN eS ET LCA TION 

The purpose of the experimental investigation was to 
(a) observe and analyze the received radial particle-velocity 
and displacement waveforms in an isovelocity layer; (b) 
compare these waveforms with theoretical waveforms produced 
as output of computer programs based on the theoretical 
developments; and (c) vary the parameter c/au) , over a 
range of values to see if the solutions broke down when 
the basic restriction c/a w , <1 was violated. The greater 
part of the experimental investigation was made using the 
step-function input; the gated-sine wave input was used 
only infrequently. Experiments were conducted using the 
procedure described in the following paragraphs. Five basic 
variables characterized each experiment; the input function 
(gated siné=wave or step-function), the range r, the 
transmitter radius a, the layer depth d, and the speed of 
sound c. 
Baepemiincweal Procedure 

Since c/aW , was a fundamental parameter of the 
solutions, the input function and the value of c/aW , to 
be used were first chosen for each run. ‘The @.5 cm 
transmitter was used for values of c/aw m erom 0.35¢Ee 
1.00; the 1.0 cm transmitter was used for values ranging 
from 1.00 to 3.00. The temperature of the water was 


measured to determine the speed of sound c in the layer (18) , 
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With c, and the desired value of c/a Ww 7 determined 
the necessary value of @, could be calevia ede tank 
was then filled to a convenient depth and Soeed as a 
pressure-release cavity of known length and width. ‘The 
depth of water was adjusted until the maximum resonance 
peak of ne fundamental mode was observed at the 
desired &,. From WG, and the known length and width of 
the tank the effective depth was calculated for the computer 
Re oer atii which was written with depth rather than Wy, as 
an input. ‘The end plate and transmitter were placed in 
the tank, oi input signal was fed to the transmitter, and 
aire este was placed at the range of 61.5 cm (used in 
all experiments) by measuring the time-of-flight of the 
signal from transmitter to receiver. 

The received waveforms were photographed from the 
Mee iloscope display. For the step input, several photo-= 
graphs were taken of each waveform. An expanded time 
scale was used for photographs of the leading portion of 
waveform where the high frequencies were contained; a 
more compressed time scale was used for sSuceeding photo- 
graphs. For a typical experiment, three photographs of 
each waveform were made utilizing, for example, 5 ~rsec/ 
em, 10 ~rssec/em, and 50 s~ssec/cm time scales. For the 
step input experiments, both received particle displacement 


and velocity were photographed in most cases. The photo- 


oe 


graphs were then read using dividers and a scale to determine 


the values of delayed time oe at which axis crossings 
occured. ‘The amplitudes of the peaks in the waveform were 
measured in a Similiar fashion. Sample photographs are 


ae in Big. 4.1. Data read from all photographs are 
Seneanend linea opencix i. 

For each experimental run, a corresponding computer 
run was made using the same values for r, a, c, and d 
(and hence c/a =) as inputs to the SacBee pregram. 
As has been shown in Section 2, a solution for particle 
pees ose or velocity in the layer can be expressed as 
the sum of the corresponding duct solution and a "layer 
Bunwectm ‘The computer programs could rierenee be 
written so that their output for each run contained a 
series of graphs of the duct solution, layer solution, and 
"layer correction" waveforms. The program for the step 
input produced a series of graphs for both particle dis— 
placement and velocity waveforms. The program for the 
gated=sine input.graphed-the vélocitphseries only. Comparison 
of graphs within a series provided convenient visual evidence 
of the relative importance of the "layer correction" with 
respect to the duct solutions as a function of the parameter 
c/a ud n° “4+ portion of the step input program was used 
to find the axis crossings of the theoretical displacement 
waveforms for the duct and layer solutions. 
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Another portion of the program evaluated the magnitude 
and delayed time of the amplitude peaks in the layer—solution 
displacement waveform. In this manner the theoretical values 
One eae: and (ee » and the magnitudes of times of the peak 
amplitudes could then be compared with the observed values 
ene a and amplitude peaks found in the experiment with 


the corresponding value of c/a Ww z. 
Analysis of Data 

Pe aie of the response of the layer to a step-= 
function input proved to be more involved than was originally 
estimated. For this reason, and because of the limited 
time espaol: for this investigation, data from the runs 
performed with a gated sine-jwave input were not analyzed. 
The following analysis pertains to the step-function inp 
data only. 

Recall that the simple theory in Section 2 based on 
the trigometric approximations for the Bessel functions 
suggested that if values of Ce are plotted versus (a. 5 a 
= eeene line should result. This relationship would offer 
advantages in data reduction because a family of axis 
crossings could be summarized by the slope of a straight 
line. For this reason, graphs which plotted the observed 
and computed values of Ce versus the computed values of 

Ges were made for each run. Figures 4.2 —- .4 are 


representative examples of these graphs. ‘The solid line of 
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unit slope represents the computed values of C. eee 


my 
circles are the computed values of Gm » and Ehescroccec ears 
the observed values of Com ». it is clear that each set of 


points can be joined by a straight line which quantitatively 
describes by its slope the predicted and observed shift in time 
of the layer solution axis crossings as compared with those of 
Eme duct. 


A graph was then made which plotted the predicted and 


observed values of oo versus c/aw,, (Fig. 4.5). The simple 
om 
theory would suggest plotting lan versus a/r (c/aw)*, how=— 
om 


ever, r was fixed for all experiments and only two values of 
a were used. Furthermore it is the validity of the layer 
solutions for various values of C/ ate which is the topic of in- 
vestigation here, rather than the validity of Eq. 2.53. 
Analysis of the amplitude peak data was also made for 
step input solutions. The experimental amplitude peaks read 
from each waveform photograph were first normalized to the 
value of the highest peak in the waveform. The predicted 
peaks were normalized in the course of the computer program. 
The ratio of these normalized quantities was then plotted as 
a function of the delayed time of the predicted amplitude 
peaks for each run. Exact agreement in the shape of the 
transient waveform between experiment and theory would 
result in a line of zero slope. Figure 4.6 shows plots for 


three values of c/awWn.- 
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Figure 4.6 


Babngetone: Analysis 


An analysis of the assumption that the semi-cylindrical 
source and steel backplate combination could be used in lieu 
of a cylindrical source was necessary. ‘The complex 
impedance of the plate was evaluated from relations in! 19) nor 

the range of frequencies contained in the received wave-~ 
forms. The impedance mis-match was sufficiently large to 
justify assuming that the plate was rigid. re directivity 
of each transmitter when mounted on the steel plate was 
investigated in order to compare the resulting radiation 
pattern with the uniform radiation assumed in the theoretical 
deveopment. Directivity patterns are Showa in Bigs. aay 
and 4.8. In order to determine the effect of the deviation 
of the actual patterns from the theoretical patterns, each 
transmitter was clamped with foam rubber so that only a 
90 degree arc in the center of the transmitter was radiating 
sound. Except for small changes in amplitude, no 
differences were observed in the received waveforms as 
compared to those received for radiation over the full 180 
degrees of the transmitter. It was therefore assumed 
that radiation from the transmitters could be assumed 
uniform and that the transmitter-plate combination could 


be used to approximate a cylindrical source. 


on. 
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The rise time cf the experimental system was determined 
to be about 0.5 psec. Therefore, instead of being a perfect 
step-function input, the velccity at the source was initially 
Heaenibed by 


on 
VG) = j-e 


where We = 2 x 10° see is the reciprocal of the rise 
time. 
The group velocity of the waveform is given by 


C= =p = ali-(@sfo)y, (4.1) 


Under the condition Ww?77a, , and with = ase y; Eq. (al 


can be written as 


= as /- £ (Prft3)* 


By rearranging Eq. 4.1, the value of delayed time at which 
a frequency a) will appear in the received waveform can be 
expressed as 
f “T [Gy & 

2 os + (2) : 
It would be expected that the amplitudes of frequencies 
greater than c€ contained in the initial portion of the 
received waveforms would be depleted up to a time given 
b 
y ae 

Cc = t(=> ) (ag 
and that amplitude distortion of the waveforms would be 
significant up to this time. However, substitution of the 
observed of and the largest value ofW, experimentally 


encountered into Eq. 4.2, shows that distortion would be 


6h, 


significant at the most for only 1.5 msec. ‘Therefore the 
rise time of the system should have no great effect on 
the received waveforms. 

It was necessary to investigate the characteristics of 
the Baa emicieks to see if the velocity of the face 
remained proportional to the step input voltage over the 
interval of time that the received waveforms were studied. 
For times much greater than the rise time, the mass of the 
diaphram can be ignored provided its velocity remains 
reasonably constant over the time interval of interest. The 
resulting equivalent circuit of the transmitter with a source 


(17) 


resistance added is shown in Fig. 4.9. Writing loop 


equations for the equivalent circuit gives 


- 9 tl icdt - 
e-ik,- af tee diz 


ee. feve dt + BVe (<<) — ¢ [4dt=o 


which enables us to write the Laplace transform of the 


diaphram velocity as 


ae se 
ae a ® 
Vv. (Sk: A 7 s* + (SEs 4 Es s +(¢ am (4.3) 
Pr Ce eS, $ Lo 
GC is determined from the electro=mechanical relationships 


ei?) 


of the transmitter and is defined as 


Ce E, 


ge ee 
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Nominal equivalent circuit vaiues for the i cm radius 
transmitter with a <0 volt step input applied are shown 
in Pig. 4.9. Substituting these values into Eq. 4.3 


and Weng the inverse transform yields 
-4 —~340t 
Vi. = 34 xo e 

Since the mass of the transmitter diaphram is ignored in 
this approximation, the particle velocity is a maximum 
eae after the application of the step input and then 
decays exponentially. Substituting t=360 Aisec (the longest in- 
terval ofdelayed time over which received waveforms were 
observed) into Eq. 4.3 shows that the velocity has 
ee ecied Coma percent of its maximum value at the end 
of this time. The tailing off of the amplitude peaks of 
the Keeemee waveforms, as shown in Fig. }.6, is consis tema 
with the decay in velocity of the transmitter diaphram. 

UnGertaimeles Inna. ‘c, r,Q),, and (han were the 
possible sources of error in this investigation. The effective 
radius of each transmitter was measured to the nearest 
1/2 mm, which could have introduced an uncertainty of + 2 
percent in the value of a for the 2.5 cm radius transmitter, 
and + 5 perceme ime thesvaluc of@a ior che 1.0 cm radius 
transmitter. Measurement of the temperature of water 
to the nearest degree restricted the (systematic) error in 
Ceeemeemichin + 0.5 percenes “The methed@ot determining 1; 


the time-of-flight of the signal from transmitter to 
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ce ohm 


Rs 


Cy = 1000 pf 


A= 4.8 x 1077 me 


fre =1.5 x 10° ke/sec-me 
Cg = 107° m/N 
Bg = 300 volts 


Xp = 3X 107? m 
p = 10-¢ N/volt 
e(t) = cO volt step input 


Equivalent circuit of transmitter 
Figure 4.9 
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fm eolvers fave 4a possible systematic error of + 0.5 percent, 
which in turn resulted in an estimated error in r of eee EES: 
percent. 

‘Two sources of error were possible in the determina- 
tion ofw a by excitation of the tank as a pressure-release 
easieos Speewwedmemrerularities of + V/iemm in thes pottom 
surface of the tank could introduce an error in the local 
effective depth of the layer. Since experiments were 
always conducted in the same portion of the tank, a 
systematic error in Ww would result. ‘This error could 
range from + 0.3 percent’ (for the 9.42 cm depth) to + 1 
pencent (for the 2.80 cm depth). The random error 
possible in determining the frequency of the resonance peaks 
of the standing wave spectrum proved to be insignificant 
relative to the possible systematic error. 


The experimental values of (. could be influenced by 


mv) 
Sy secMocleweuuows ON@ehe Onder OF tmleeecr cont caused jan, 
nonlinearity and miscalibration of the oscilloscope time-base. 
As mentioned previously in this section, axis crossing times 
were read from a series of photographs with different 
time scales. This method was used so that the initial 
portion of the received signal would always be at the left- 
hand edge of the picture and the axis crossings of interest 
would generally occur in the right half plane of the display. 


Thicepaeccdure resulted in a random errompot + 0.5 percent 
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on all time scales since the display could be read to within 
1/h scale division out of 50 divisions. 
he errors in the computed values of ig. were influenced 

predominately by the error insect ice mace n 2ppears 
in the argument of the Bessel functions { Ww, V “-T* ) from 
which the ‘me 's are computed. The values of lager, were 
determined when the argument equaled any one of a set of 
constants Ge which resulted in a zero value for the 
amplitude of the particle-displacement solution. ‘Thus the 
argument can be written as 

w Veer = Lo, We, -T\ 4,47) = wW,, V tn Gaon = Cm 
and if we have (ae eK Ay » then 

Var t,2% = Cm 


1 


Therefore the + 1 percent error in W dé could lead to + 2 


and 


percent error in the computed values of TT, 
The errors in the slopes of the curves fitted to the 
eeSinental e's in Figs. 4.2 = .4 were determined mainly 
by the systematic errors in these values of 7, because 
the random errors were minimized by fitting a straight line 
through the points. Therefore the error in each slope 
boas CSCimated to be + 1 percent. ~@here was a systematic 
error of mee Percent) nett eucompMeccmalves Ot gee > hence 
the error in the slope of the computed lee versus 


om 


lines was + 2 percent. 


69 


The errors in c, a, and w, being independent can be 
ee aniied by the root-mean-square law to yield the error 
ic) 3) Glee) Ee penecni error in c/auy , varied according 
to the values of a and W, involved. Table 4.1 shows the 
result of this combination for representative values of 
c/a Ww: nae flags on Figs. 4.2 - 4.5 indicate the 
relative magnitudes of errors in the experimental and 
computed values of Ge and slopes, and the relative 


magnitude of error in c/a W i 


Table 4.1 


‘Tabulated Prrors for Quantities Displayed on Figures ).2-4.5 


Percent Prror 





70 


oe RESULTS AND CG@nN Geeta N] 

The investigation has shown that for c/a Oo, <1, 
the predicted and experimentally observed shift in axis 
crossings agreed to within + 2 percent. Where the 
restriction c/aw < 1 was violated, the theory failed to 
predict the observed shift in axis crossings. As c/aw, 
was increased the failure of the theory became more 
pronounced. 

The correlation between the predicted and observed 
amplitudes was less favorable. The divergence between 
the predicted and observed amplitudes in the initial portion 
of the received signal was unexplained. ‘The decrease in 
the observed amplitudes in the latter portion of each 
Peesived signal is not inconsistent with the roll-off of the 
step-function input in velocity. 

Therefore, with the exception of amplitudes, we can 
say that the theoretical solutions correctly described the 
propagation of simple transients in an isovelocity layer with 
perfectly-reflecting boundaries as long as the parameter 
c/a Ww was &1; and that, although the theory failed for 
values of c/aW , 7 1, it did correctly predict the trend 


that the observed waveforms followed. 
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Pies IN IT X A 


A-1: Summary of the Propagation of Transients in Waveguides. 
Results for the propagation of simple transients (step- 
function and gated sine-wave input) in ducts which have 


(4,10) 


appeared in the literature and which will be of value 
in this research are presented below along with some addi- 


tional Laplace transform pairs that are needed. 





jay ye byee gaa A 


Assume a waveguide of uniform cross-section, infinite 
length, and perfectly-reflecting boundaries, which is filled 
with a lossless fluid and excited by a velocity souce at z = 0. 
See Fig. A.1 The wave equation for the acoustic velocity 


potential D can be written in cylindrical coordinates as 


ov, Ld >> 1 > 
(Se rer 557 ot 2.) $ (143,t) =O (A.1) 


an 


where r is the radial direction, z is the longitudinal direction, 
and radially symmetric processes have been assumed. 


The separation of variables 


D(v3,t) = 2, D,, (v,33¢) = di Racr) B63, t) (A.2) 


leads to the solution R= J (@ r/c) where J, represents 
n 
the zeroeth order Bessel function of the first kind, and 


the partial differential equation 


s  _/o,\*_ 1 O° A.3) 
ey (Zh Se, | ee) —— 

which describes the propagation in the 2 direction. 

The value of the Separation constant eae depends on the 


nature of the boundary at r = a (the wall of the duct) and 


the cross-sectional area of the tube. Because the transverse 


properties of the duct affect only Ww and not the form of 
the solution of Eq. A.3, all functions of r have been 
suppressed. They can be reintroduced with no difficulty 


when desired. 


The Laplace transform of Eq. A.3 with quiescent initial 


conditions is 





2 ey =. oe 
2 Z,= Sn Z (Av) 
O23. c* é 


= 


where the bar over the function indicates the Laplace 
transform of that function, and D, _ is a fragment of the 
velocity potential containing z and t. Solving Eq. A.4 with 


the boundary condition 8, = B (0,t) gives 


B= G exe lr or] as) 


where T = z/c is the time of flight of the signal from the 
source to the receiver. 
The relation 4.,, =f, aPn is used to obtain the 
03 
Laplace transform of the fragment of the longitudinal 


particle velocity: 


a. UL me / se (ane 
ata = “7p (r) = <2 S*4Hn” exp LT Hor | ° \ ae 
é . n ral 





Application of the boundary condition on velocity at the 


sKoneaetste) vie (01) Vite) casts Eq. A.6 into the following form: 
Van = V. exe -TVSHOr | (Av 


The relation 





is evident from examination of Eqs. A.6 and A.7. If 
this value of B. is substituted into Eq. A.5, the Laplace 


transform of 2. in terms of the velocity input becomes 


ry a ~eV. exe[-TVs%ar| 
EO 


The fragment of the longitudinal acoustic pressure is 


(A.9) 


obtained from the velocity potential by means of the 


i = a i - = ) e i | 
expression Pa n= Rr) ‘4 r This corresponds to 
a = - a sf in the Laplace transform domain. We 


therefore have the result 


= ecsV,_ exe(-TVs%0,] (A 40) 


Vs a 


The fragment of longitudinal particle displacement may 
be obtained either by integrating the fragment of longitudinal 
particle velocity in the time domain or by multiplying by 1/s 
in the Laplace transform domain. With the latter alternative, 


the particle displacement fragment S4n can be written as 


6 = Vi exe LT VS905] (A.11) 
nm |UOUS 


77 


Fiquations A.7 and A.9 = A.11 are the Laplace transforms of 





the acoustic variables of interest. The inverse transforms 
for a step-function input ( (Vi, = 1/s) and a gated sine-wave 
= UW ; 

input ie mecneyeo ) are as follows: 


For a step-function input 


@= -2 at eye 
hn My oly i) nt (A.12) 
are 
RG ay 
J,(are) + a2 (=) J, CA86) ( Agi 
ale =¢c_), (ARC) (A.1)) 


A 


g — & -) (2m- “(8 ae, (nee) (A.15) 


mn>) 


Where ARG = W “ ie _pe 
For a gated sine=-wave input: 
~ L¢ 
?,= “PIU, (w,aee) - U. (nv-, aes) | (A.16) 
LE U, (my, RES) +p U, (av. Ree) (A.17) 
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¥3,= CCr | U, (wae) ~ U, (uw-, nee) | (A.18) 


mt 


..* -+ ilu, (i, nes) —U(rs-, ARG) -2F (EN J cree (A.19) 


Where ARG = Uden t -T 


ws (4-T)(1# Ge) 


The U's represent Lommel functions of two 


variables 


+2m 


On (nz, ARC) = a _(-!) "( :) T, 4 ee) 


A-2: Development of Additional Laplace Transform Pairs 





for Step Input in. Velocityein ame bsovelocitwweaver: 





The transform pairs developed in this section are 
those needed for the layer solution in addition to those 
summarized in section A-1. EF denotes the transform under 
consideration; f(t) is the inverse transform of F. All f(t)!'s 
are defined as zero fort £& T. 


Hh 


ic = sV. exe (- TVs," | (A.20) 


S24W we 


eg 


With a = 1/s, Eq. A.20 becomes 





F= exp -TVst0,* | (AT2m) 
ne 
Sua 
Equation A.21 is of the same form as Eq. A.7 with 
a Fe 
Vo = Fee and W= td, . Therefore the inverse 
transform of Eq. A.7 may be used with wWe=tw, . The 


inverse transform of Eq. A.7 is 


D, (as, nee) + U, (ms, ARG) (A.22) 


where ARG = wD, oe 


and = Usp = W(t-T)(/t Ye ) 


‘The quantity C q /c can be written as ae n/ud 
q / | — (Aris), 


With (Us WwW, , we have Avy= UW = w _(t-T). 
This in turn means that both terms of Eq. A.@é are 


C44) 


equal. Using the relation 


aU, [o,(4), Ae6| = 2 Le ae olen) ees (A.23) 


Sate ae vtam 


the desired inversion is obtained. With \) ari 


ARG a t©-Té, Fq. A.21 may be rewritten as 


£(e) ZU (ns, Ars) = ay ')" le Ve “J (ee) (A.2h) 
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= V, exe (-TV S*W,* ] 
F S \J sted, 


(A.25) 


With Vo = 1/s, Eq. A.25 becomes 


ie oa exp(-T (s4+o),*| (A.26) 


S*\ ste, 





The identity 





| IVv+l 
a (A.27 
Xx eK? = 
is used to generate the desired inverse transform. Squaring 


Eq. A.27 gives 


(v+1) 
(+) 4 y (V+!) Fook - Aiea 








Letting x = Wn/s in Eq. A.28 and rearranging yields 
sand I(t!) 
ra 
Sener On 
—- SS +\ 


Solving Eq. A.29 for i/s and substituting the value obtained 


into Eq. A.26 results in 


Fed {ud (ve) Foe 7" exel-TV sae} ae 
Vstewn 





oie 


8) 


The inverse transform of Eq. A.30 is aueeeres . There-— 


fore with \)Yy+t1li=m , 
£, 4 < 4-TY” 
(t) = se >, m (eT) J. (aee) (A.31) 


A-~3: Inverse Laplace Transforms for a Gated Sine-Wave 
Inputein Velocity in an Isovelocity Layer. 


The following acoustic quantities could not be 
completely solved. Although the expressions listed are 
incomplete solutions, it is worthy of note that the leading 
terms are the corresponding terms of the duct solutions 
in Section A-1. In the equations listed tt f denotes 


the inverse Laplace transform of the quantity in brackets. 


C 2 ES \ Stu? 
fe $= SCH el 00 t HAY48 af (aac cede d (A.32) 


(Strust¥st+-os,2) 


(stew*) (524w),2) 


Ee elu, (aw, res) - U, (0500) 4 +e ys nd eh (A.33) 
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Sie ace ale S /t- 
\E . 2 [Us nea) +05 (spn) +25 (5) "T (nes) (A.34) 


- Pee ae (: a)T [ene] 


where ARG = Wy, Veer? 


pss = w(t-TY1tS/e) 


U 


; Uy A 
yiw, Are) = Lommel Function of two prdics 5 oe 


Uy (ave, ARG) = 5 ney y J (ane) 
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APPENDIX 8B 
EXPERIMENTAL DATA 


DISPLACEMENT DISPLACEMENT VELOCITY 
AXIS CROSSINGS AMPLITUDES AXIS CROSSINGS 
TPRI IN USEC (NORMALIZED) TPRI IN USEC 


STEP INPUT RUN $3 
R=6125CM A=2e54CM D=2e80CM 


C/AWN=0 235 


020 
1210 
2285 
2250 
9200 
13.3 
18.5 
2463, 
3028 
38.20 
4509 
24e2 
63.8 
7365 
84.0 
9640 


STEP INPUT RUN Si2 


C/AWN=0 052 


020 

2.50 
2e10 
10.7 
lve? 
2601 
366 


48.4 
6163 


1620 
Sle? 


10865 
126.20 
14540 
164.20 
18465 


R=612e5CM 


Ool75 
-00379 
00442 
-00592 
02660 
-0e718 
0«864% 
—Oaaec 
12000 
-0-898 
02893 
-02718 
0Oe718 
-0e592 
02563 


A=2e54CM D=4e15CM 


02198 
-0 2396 
02636 
-02833 
12000 
-0.972 
02905 


-0 +860 
02838 


-0 2838 
Oe72l 


-02694 
02622 
-02613 
02532 
—00477 
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TMP =20DEGC 


095 
2035 
435 
7265 
113 
1665 
21le7 
2802 
3403 
4204 
50 «3 
59 e2 
68e2 
195 
91-20 


TMP=20DEGC 


NOT RECORDED 


85 


DISPLACEMENT DISPLACEMENT VELOCITY 
AXIS CROSSINGS AMPLITUDES AXIS CROSSINGS 
IPRie CNBUSEG (NORMALIZED) TPRI IN USEC 
STEP INPUT RUN S6 
C/AWN=0 «69 R=61e¢5CM A=1e00CM D=5e50CM TMP=20DEGC 
020 00414 220 
3450 -0.68/7 6.2/0 
9260 02990 1464 
18.4 -1-000 2302 
2948 0.919 3806 
45.8 -0.909 52¢5 
6440 0-828 1305 
85.0 —-02854% 930 
1060 02/798 119.0 
131.0 Oita 140-0 
168.0 02768 169.0 
183.0 -0+-687 193.20 
212-0 02687 
24020 —-0665/7 
26840 02667 
2980 -0-581 
32940 02616 
3600 
STEP INPUT RUN S13 
C/AWN=0 2-85 R=61-¢5CM A=2e54CM D=6e78CM TMP=21DEGC 
020 00357 NOT RECORDED 
4.20 -02697 
1326 0-871 
2603 —-02874 
44.5 00994 
6569 -0.931 
902 Oel77 
118.0 —-0-2834 
14765 1-000 
17965 —0 2943 
21265 0.783 


DISPLACEMENT DISPLACEMENT VELOCITY 
AXIS CROSSINGS AMPLITUDES AXIS CROSSINGS 
TPRI IN USEC (NORMALIZED) TPRI IN USEC 


STEP INPUT RUN 314 


C/AWN=1200 R=61e¢5CM A=2e54CM D=7e¢97CM TMP=21DEGC 
0.0 02538 NOT RECORDED 
7200 -0.2813 
19.0 0869 
3620 -0+875 
Doe 5 02943 
89e2 -02956 
12020 12000 
1560 -0.-863 
1900 0-769 
235.20 -0.2619 


STEP INPUT RUN S4 


C/AWN=1 00 R=61¢5CM A=1eO00CM D=3e14CM TMP=20DEGC 
020 Oo176 0-80 
1.25 —-0e412 2044 
3045 02559 4.90 
6640 -02598 8/5 
1065 02696 133 
1620 -02814 289 
2129 0.961 woec 
2802 -1.000 3204 
35e7 12000 405 
4502 -0.971 49e2 
D365 0922 580 
6365 -0 +863 6865 
74.0 0784 71920 
85.0 -02/706 9125 
98.5 0667 
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DISPLACEMENT 
AXIS CROSSINGS 
IN USEC 


TPRI 


DISPLACEMENT 
AMPLITUDES 
(NORMALIZED) 


STEP INPUT RUN S9 


C/AWN=1.50 


0.0 

2230 
6.50 

12.8 

2066 
30.6 
EDS) 
515 8. 
alee S 

88.0 
107. 0 
Sree 
145 0 
1690 
191.0 
2A SO 
23920 
262.20 
2860 
312.0 


R=61e¢5CM 


STEP INPUT RUN SSS 


C/AWN=2 200 


OO 
3.80 
eye 
Ziel 
3464 
52.5 
Wize 
95.5 
119.0 
145.0 
176.0 
203-0 
235-40 
2650 
29520 


R=612e5CM 


A=1e00CM 


02307 
-0-2586 
0.793 
-—0 2923 
1-000 
-02964 
0.979 
—-0-943 
0-821 
—-0-686 
0.750 
-02700 
02.5923 
-0 2543 
02521 
-02543 
02564 
—02479 
02486 


A=] e00CM 


Oe Jou 
-0 2649 
1-000 
-1.000 
0.979 
-0.918 
0.907 
-0.701 
Oi ol 
—-0-639 
0.639 
-0 639 
02639 
-0.-598 
02495 
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VELOCITY 


AXIS CROSSINGS 
IN USEC 


TPRI 


D=4-e71CM 


D=6e28CM 


TMP=21DEGC 


1.30 
420 
9220 
165 
250 
3507 
490 
650 
1940 
980 
116.0 
13210 
154-0 
178.0 
202620 
225-40 
246-0 
26920 
297.20 


TMP=20DEGC 


2210 
7260 
1565 
2602 
41.20 
5965 
1945 
101.0 
127.0 
155.0 


DISPLACEMENT 
AXIS CROSSINGS 
LPR iy wUSEC 


STEP INPUT RUN 
C/AWN=3.00 


020 
7250 
21.0 
40.5 
67-20 
100.0 
13260 
184.0 
2220 
272.20 
312.0 


DISPLACEMENT VELOCITY 
AMPLITUDES AXIS CROSSINGS 
(NORMALIZED) TPRI IN USEC 


or 
R=61e¢5CM A=1le00CM D=9442CM 


02432 
-0 «892 
12000 
-0 6826 
0822 
-0./784% 
02526 
—-02427 
02376 
ORE he 
02639 
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TMP=20DEGC 


4.50 
15-0 
3140 
5360 
84.5 
11540 
150.0 
210.0 


VEEOG lay 


AXIS CROSSINGS 
IN USEC 


TPRI 


GATED SINE INPUT RUN Gl 
C/AWN=1 200 


380 
9-60 
Ga =, 
2leo 
2643 
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